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We show how to calculate orbit implications (based on 
Nielsen- Thurston theory) for horseshoe type maps arising in 
chaotic time series data. This analysis is applied to data from 
the Belousov-Zhabotinskii reaction and allows us to (i) predict 
the existence of orbits of arbitrarily high periods from a finite 
amount of time series data, (ii) calculate a lower bound to the 
topological entropy, and (iii) establish a "topological model" 
of a system directly from an experimental time series. 
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I. INTRODUCTION 

Braids arise as periodic orbits in dynamical systems 
modeled by three-dimensional flows (see Fig. 1) |^|-Q. 
The existence of a single periodic orbit of a dynamical 
system can imply the coexistence of many other peri- 
odic orbits. The most well-known example of this phe- 
nomenon occurs in the field of one-dimensional dynamics 
and is described by Sarkovskii's Theorem g]. Less well- 
known is the fact that analogous results hold for two- 
dimensional systems B . In one-dimensional dynamics it 
is useful to study the period (or the permutation) of an 
orbit [^ . In two-dimensional systems it is useful to study 
the braid type of an orbit 0] . Given this specification, we 
can ask whether or not the existence of a given braid 
(periodic orbit) forces the existence of another; as in the 
one-dimensional case, algorithms have recently been de- 
veloped for answering this question [0-^ . 

As originally observed by Auerbach et. al., unstable 
periodic orbits are available in abundance from a single 
chaotic time series using the method of close recurrence 
[^,Q. By a "braid analysis" we propose to analyze 
a chaotic time series by first extracting an (incomplete) 
spectrum of periodic orbits, and second ordering the ex- 
tracted orbits according to their orbit forcing relation- 
ship. As shown in this paper, it is often possible to find 
a single periodic orbit (or a small collection of orbits) 
which forces many orbits in the observed spectrum. This 
single orbit also forces additional orbits of arbitrarily high 
period. This analysis is restricted to "low-dimensional" 
flows (roughly, flows which can be modeled by systems 
with one unstable Lyapunov exponent), however it has a 
strong predictive capability. 

We would also like to point out that this analysis gives 
us an effective and mathematically well deflned "pruning 
procedure" for chaotic two-dimensional diffeomorphisms 



[]l2| . Instead of asking for rules describing which orbits 
are missing (pruned), we instead look for those orbits 
which must be present. For low period orbits (say up 
to period 10) this procedure can predict all those or- 
bits which must be present in the flow. This procedure 
will usually miss orbits of higher period, however from 
an experimental viewpoint the low period orbits are the 
most important and accessible. Orbits of low period of- 
ten force an infinity of other orbits. This is illustrated in 
one-dimensional dynamics by the famous statement "pe- 
riod three implies chaos" U^. An analogous statement 
in two-dimensional dynamics is that a non-well-ordered 
period three braid implies chaos jlj] . 

This paper is organized as follows. In section II we 
briefly review the theory of orbit forcing in one- and two- 
dimensional maps (three-dimensional flows). In section 
III we show how the braid analysis method works by ap- 
plying it to times series data generated from the Rossler 
equations. This section also discusses some useful nu- 
merical and symbolic refinements to the method of close 
recurrence. In section IV we apply the braid analysis 
method to data from the Belousov-Zhabotinskii reaction. 
Our analysis builds directly on the original topological 
analysis of this data set due to Mindlin et. al. |^5[. In 
section V we offer some concluding remarks by indicating 
how this method can be extended to low-dissipation and 
conservative systems. In the examples studied in this pa- 
per we do have good control of the symbolics. In princi- 
ple, though, this method does not require good control of 
the symbolics (a good partition) and can thus overcome 
some of the current difficulties associated with finding 
good symbolic descriptions for (nonhyperbolic) strange 
attractors |16| . 



II. ORBIT FORCING THEORY 

In this paper we study in detail the orbit structure of 
horseshoe type attractors. In particular, we will make 
a good deal of use of those features of one-dimensional 
dynamics which must always carry over to two dimen- 
sions p^ ]. Therefore, let us begin by summarizing some 
well-established [p 18 1 definitions and results pertaining 
to the periodic orbit structure of maps of the line. Let 
_F : R — > R be a continuous map, and R = {xi, . . . ,Xn} 
be a period n orbit of F with its points labeled in such 
a way that xi < X2 < ■ ■ ■ < Xn- Then there is a cyclic 
permutation vr = vr(i?, F) associated to R, deflned by 
F{xi) = x-^^i). A relation <i on the set C of all cyclic 



permutations can be constructed as follows: if cr,T € C, 
then r <i (J if and only if every map _F : R — > R which 
has a periodic orbit R with tt{R, F) — a also has a peri- 
odic orbit S with 7r(S', F) = t. If this is the case then we 
shall say that the permutation a forces the permutation 
T. The relation <i is a partial order which can be cal- 
culated, in the sense that there is a simple algorithm to 
determine whether or not r <i cr for any two cyclic per- 
mutations T and a: unfortunately this algorithm takes a 
long time to execute for orbits of high period, and it is 
difhcult to use it to analyze the global structure of the 
partially ordered set (C, <i) Q. There is, however, an 
interesting subset of C onto which the restriction of <i 
is well-understood. Let UM C C be the set of unimodal 
permutations: that is, the set of permutations of periodic 
orbits of the tent map (or one-dimensional horseshoe) 

3a; < x < 1/2 

F{x) = { 3(1 -x) l/2< a; < 1 
otherwise 

(equivalently, UM is the set of permutations which can be 
realized by periodic orbits of unimodal maps whose turn- 
ing point is a maximum). The partial order <i restricts 
to a linear order on UM which is described by kneading 
theory p. 

This is a statement of topological universality in the 
one-dimensional context: there is a universal order in 
which the periodic orbits are built up in families of uni- 
modal maps of the line. In Appendix A we provide a 
listing of the unimodal order for orbits with periods 1 
through 10. It is important to understand that the state- 
ment '(UM, <) is a linearly ordered set' gives constraints 
on the possible periodic orbit structure of all maps of 
the line, not just unimodal ones: similarly, the results 
which we shall present here are applicable to arbitrary 
orientation-preserving homeomorphisms of the plane, not 
just to those which are horseshoe-like. 

A periodic orbit i? of a homeomorphism of the plane 
can be described by its braid type |Q]; the reader who 
lacks the mathematical background necessary to appre- 
ciate the following definitions should simply regard the 
braid type as a two-dimensional analog of the permuta- 
tion, and the term 'pseudo-Anosov' as one which, while 
essential for the accurate statement of theorem 1, does 
not appear in later statements. Very roughly, braids 
can be divided into three types, the so-called finite or- 
der braids — braids which only force the existence of a 
finite number of other braids (periodic orbits), — pseudo- 
Anosov braids — braids which must force the existence 
of an infinite number of other braids, and reducible braids 
— braids which can be decomposed into distinct compo- 
nents which fall into one of the first two types. Pseudo- 
Anosov braids play an important role in the study of 
chaotic three-dimensional fiows since, at the topological 
level, they force the coexistence of an infinite number 
of other periodic orbits. The term pseudo-Anosov arises 
within the context of the Nielscn-Thurston classification 



of isotopy classes of surface automorphisms, which is the 
main tool used in the proofs of the results presented here: 
the reader interested in this powerful theory could consult 
Ref. [gO|, in which there is also an extensive bibliogra- 
phy. An exposition of the relevance of Nielsen-Thurston 
theory to two-dimensional dynamics, with particular ref- 
erence to braid types, can be found in Ref. B- 

Let / and g be orientation-preserving homeomor- 
phisms of the plane which have periodic orbits R and S 
respectively. We say that {R, f) and {S, g) have the same 
braid type if there is an orientation-preserving homeomor- 
phism h ■.'B? \R ^ Yl^\S such that h o /|r2\^ o h^^ 
is isotopic to g|R2\5. Having the same braid type is an 
equivalence relation on the set of pairs (i?, /), and we 
write bt(i?, /) for the equivalence class containing (i?, /), 
the braid type of the periodic orbit R oi f. We say that 
the braid type bt(i?, /) is of pseudo-Anosov, reducible, or 
finite order type, according to the Nielscn-Thurston type 
of the isotopy class of/ : {S^,R\J{oo}) -^ {S^,R\J{oo}) 
(see Ref. |2^| for a statement of Thurston's classification 
theorem for homeomorphisms relative to a finite invari- 
ant set). By analogy with the one-dimensional theory, 
we define a relation <2 on the set BT of all braid types 
as follows: if /3, 7 G BT, then 7 <2 /3 if and only if every 
homeomorphism of the plane which has a periodic orbit 
of braid type (3 also has a periodic orbit of braid type 
7. If this is the case, we shall say that the braid type 
/3 forces the braid type 7. The relation <2 is a partial 
order p3]: as in the one-dimensional case, there is an 
algorithm for calculating this order (the so-called "train 
track algorithm") Q, but it is difficult to use it to ana- 
lyze the global structure of the partially ordered set (BT, 
<2). Still, using a version of this algorithm due to Bestv- 
ina and Handel 0|, it is possible to calculate the orbits 
forced by an individual braid up to say period 10 or so 
by hand, and we have done this for several low period 
pseudo-Anosov horseshoes braids in order to find their 
associated topological entropies (see Table |). Also, an 
automated version of this algorithm promises to extend 
these calculations to higher period orbits [ p3[ . 

We restrict our attention to the subset HS C BT of 
horseshoe braid types: that is, those which are realized 
by periodic orbits of the horseshoe map / |2J] . The de- 
termination of constraints on the order in which periodic 
orbits can be built up in the creation of horseshoes is 
tantamount to the analysis of the partially ordered set 
(HS, <2); and an important observation in discussing 
the topological universality of two-dimensional systems 
is that, unlike (UM, <i), the partially ordered set (HS, 
<2) is not linearly ordered. 

The connection between the ordered sets (UM, <i) 
and (HS, <2) is that periodic orbits of both the tent 
map and the horseshoe have a symbolic description. We 
associate to points in the non- wandering set of each map 
an infinite sequence of O's and I's in the usual manner, 
and describe a period n orbit R by its code cr, given 
by the first n symbols of the sequence associated to the 
rightmost point of R. We shall identify periodic orbits of 



the two maps which have the same code, using the same 
symbol R to denote each orbit. Such an identified pair R 
will be referred to as a horseshoe orbit (or simply as an 
orbit): the fact that the orbit is periodic will always be 
assumed. We shall write c^ for the semi-infinite sequence 
crCrCr . . . (which is associated to the rightmost point of 
R, regarded as a periodic orbit of the tent map); and 
°°c^ for the bi-infinite sequence . . . cjicrCr . . . (which is 
associated to the rightmost point R, regarded as a peri- 
odic orbit of the horseshoe). 

It will be convenient for us to make another identifica- 
tion. If i? is a period n orbit, then denote by cji the code 
obtained by changing the last symbol of cr. If cr is not 
the two-fold repetition of a sequence of length n/2, then 
it is the code of a period n orbit R' which can be shown 
to have the same permutation and braid type as R. We 
identify R and R' in what follows, using the same sym- 
bol R to denote both, referring to them as a single orbit. 
The code cr of the identified pair will always be chosen to 
end with a zero (this is to ensure that the sequence c^ 
contains the group of symbols . . . 010 . . ., which is nec- 
essary for the algorithms which follow). After making 
this identification, distinct horseshoe orbits have distinct 
permutations: however, it is possible for several different 
orbits to have the same braid type. With this in mind, 
we say that a function defined on the set of all horseshoe 
orbits is a braid type invariant if it takes the same value 
on any two orbits of the same braid type. 

Given two horseshoe orbits R and S, we shall write 
S <i Rif and only if 7r(S', F) <i 7r(i?, F); and we shaU 
write 5 <2 i? if and only if bt(S', /) <2 bt(i?, /). We re- 
fer to <i and <2 as the one- and two-dimensional forcing 
orders respectively. Our aim is to analyze <2 using the 
well- understood properties of <i. Notice that the orders 
<i and <2 cannot strictly be regarded as being defined 
on the same set: if R and S are two horseshoe orbits 
which have the same braid type but different permuta- 
tions, then they must be regarded as being equal when 
considering the order <2, but as being distinct when con- 
sidering the orbit <i. For example, the orbits /3X2 and 
Sg in Table || are two such orbits (for an example analo- 
gous example in the orientation-reversing Henon map see 
Ref. H). 

We say that a horseshoe orbit R is quasi-one- 
dimensional (qod) if for all orbits S we have R >i S :=> 
R >2 S. As proved in Ref. [0, the following result 
provides a simple classification of the qod orbits which 
have pseudo-Anosov braid type. Let Q denote the set 
Q n (0, 1/2) of rationals lying (strictly) between and 
1/2. Given such a rational q — m/n (in lowest terms), 
we write Pq for the period n-\-2 orbit which has code Cq — 
lO'^iii)i2Q>^2(q)i2 ^ ^ ^ i2o«^™(9)io^ where Ki{q) = [l/q\ - 1, 
and Ki{q) = [i/q\ — [{i — 1)/?] — 2 for 2 < i < m (here 
[x\ denotes the greatest integer which does not exceed 
x). Thus, for example, C2/5 = 1011010. Then we have: 

Theorem 1 q < — > Pq is a one-to-one correspondence be- 



tween Q and the set of quasi-one-dimensional horseshoe 
orbits which have pseudo-Anosov braid type. Moreover 
Pq >2 Pq' ■^=>q'>q for aU q, q' € Q. 

It is a consequence of quasi-one-dimensionality that 
Pq >2 Pq' 4=> Pq >i Pq' , and therefore the braid types 
of the orbits Pq form a linearly ordered subset of BT. The 
final statement of the theorem tells us that the ordering 
within this subset is simply the reverse of the usual or- 
dering on Q. 

Theorem 1 allows us to quickly identify a very useful 
subset of pseudo-Anosov braids. Some of the low period 
qod (pseudo-Anosov) orbits are listed in Table |l[ and 
a program is easily written to quickly generate all the 
horseshoe qod orbits of any desired period. Moreover, 
Theorem 1 can be used to define an invariant of horse- 
shoe braid type: the height q{R) of a horseshoe orbit R 
is the unique element of [0, 1/2] with the property that 
for all g G Q we have q < q{R) =^ Pq >i R and q > 
q{R) => Pq^iR (thus qiR) = sup{q e Q : Pq >i R}). 
Because the orbits Pq are qod, it follows immediately 
that q < q{R) =^ Pq >2 R- In fact, it can be shown 
that the second property in the definition of height also 
holds for the two-dimensional forcing order: that is. 



Theorem 2 Let i? be a horseshoe orbit and q G 
Then q < q{R) =^ Pq >2 R and q > q{R) = 
In particular, height is a braid type invariant 



Pq^ 



Q 

2R- 



We can also easily find the finite order braids because 
the finite order braids of the horseshoe are exactly 
those whose period equals the denominator of the height 
|17| , p6[ . Because the height is defined in terms of the 
one-dimensional order, it can easily be calculated using 
kneading theory. We first define the height q{c) G (0, 1/2] 
of a semi-infinite sequence which begins c — 10 ... , and 
which contains the group of symbols . . . 010 .... To do 
this, write c in the form c — lO'^^l^^O'^^l'^^ . . ., where 
Ki > and iJi is either 1 or 2 for each i, with /i; = 1 only 
if Ki+i > (thus the Ki and ^i are uniquely determined 
by c, and the fact that c contains the group . . . 010 . . . 
means that fis — I for some s). Let 



Iric) 



2r- + El=i«.'(2r-l) + E[=i^ 



for each r > 1, and let s be the least positive inte- 
ger such that either fis — I or r\^^lli{c) = 0. Then 
n|=i/i(c) = (x, y] for some x and y: we define q{c) = y if 
/is = 1 or Is+i{c) > r\f^i{c), and g(c) = x if /is = 2 and 

/s+i(c) < n,Li/^(c). 

Theorem 3 q{R) = <i{cr)'- in particular, q{R) is posi- 
tive and rational. 

For example, let R be the period 20 orbit with code 
Cfl ^ 10000110001100001100. We have ki = 4, /ii = 2, 
K2 = 3, /i2 = 2, K3 = 4, /i3 = 2, K4 = 2, and /i4 = 1. Thus 



h = (1/6,1/5], h = (2/11,2/10], /3 - (3/17,3/16], 
and h = (4/21,4/20]. Since 4/21 > 3/16 we have 
h n h n I3 n h = O, and h > /i n /2 n /s. Therefore 
q{R) = max(/i n /a fl I3) = 3/16. 

There is a relationship between height and a weU- 
estabhshed braid type invariant: the height of a horse- 
shoe is equal to the left hand endpoint of its rotation 
interval. A proof of this is given in Ref. jlj, where a 
practical algorithm for determining the rotation inter- 
val of a horseshoe is described. It can also be shown 
using the algorithm for determining the height that for 
each q = m/n (in lowest terms), Pq is the only period 
n + 2 orbit of height q; thus it is the only orbit of its 
braid type. This observation enables us to determine 
exactly which of the orbits Pq are forced by a given or- 
bit R in the two-dimensional order: we now present an 
algorithm for this purpose. If i? is a horseshoe orbit 
then we define the depth r{R) G (0, 1/2] n Q of i? as 
follows: consider all groups of the form ... OHIO .. . or 
. . . 01010 ... in the sequence °°c'^] suppose that there are 
I such groups gi, . . . ,gi contained in one period of °°c^. 
li I = then r{R) = 1/2. Otherwise, for each i < / let /j 
be the code obtained by starting at the last 1 in gi and 
moving forwards through °°c^, and bi be that obtained 
by starting at the first 1 and moving backwards. Then 
riR) = mini<i<imax{q{f^),qibi)). 

Theorem 4 r{R) is the unique element of [0, 1/2] with 
the property that for all <? £ Q we have r{R) < q =^ 
R >2 Pq and r{R) > q =» R^2Pq- In particular, depth 
is a braid type invariant. 

The following corollaries from theorems 2 and 4 will 
be very useful in our braid analysis since they allow 
to us to apply one-dimensional theory to calculate two- 
dimensional orbit forcings. 

Corollary 1 Let R and S be horseshoe orbits. If 
r{R) < q{S) then R >2 S. On the other hand, if 
qiR) < q{S) and r{S) < r{R) then R^2S and S^2R- 
thus orbits of the braids types R and S can exist inde- 
pendently of each other. 

In fact, a much stronger result can be proved: if r{R) < 
q{S) then every homeomorphism of the plane which has 
a periodic orbit of braid type R has at least as many peri- 
odic orbits of the braid type S as does the horseshoe (the 
corollary only says that every such homeomorphism has 
at least one such orbit). Thus all of the periodic orbits 
of the braid type S must be created before any of the 
periodic orbits of braid type R in any family of home- 
omorphisms leading to the creation of a horseshoe. We 
can use the following corollary to locate some of these 
orbits. 

Corollary 2 If a homeomorphism of the disc has a qod 
orbit Pq, and R is another orbit which is forced by Pq in 



the one-dimensional order (which means that its height 
is > g), then the homeomorphism must have at least one 
orbit of the braid type of R. If in addition the height 
of R is strictly greater than q, then the homeomorphism 
must have at least as many orbits of the braid type of R 
as does the horseshoe |p7| . 

For example, if a period 7 orbit with the braid type of 
1011010 (height 2/5) is extracted, then there must be 
at least one orbit of the braid type of 10110 (since this 
braid type has height 2/5), and at least as many of the 
braid type of 101110 as in the horseshoe (i.e., 2) need 
exist (since the height of these orbits is also 2/5). On 
the other hand, both of the orbits in the period 6 pair 
101115 must exist since the height of these orbits are 1/2. 
A listing of the low period qod orbits and some of the low 



period orbits which they force is found in Table III 

Since the qod orbits essentially inherit the one- 
dimensional forcing order we can use one-dimensional 
methods to calculate their topological entropies. This 
information can be used to obtain lower bounds for 
the topological entropy of a partially-formed horseshoe. 
More precisely, given a horseshoe orbit R, let h{R) de- 
note the smallest possible entropy of a homeomorphism 
having an orbit of the same braid type as R. Now we 
know that for all q > r{R) we have R >2 Pq, and hence 
the topological entropy h{R) > h{Pq). However, since Pq 
is a qod orbit, it can be shown that h{Pq) is equal to the 
entropy of Pq regarded as an orbit in one dimension: this 
can be calculated using standard transition matrix tech- 
niques. A Mathematica program based on the method of 
Block et. al. [g8| has been written which allows us to cal- 
culate the topological entropies of the qod orbits. Some 
of the results are listed in Table fU A graph of the func- 
tion {q,h[Pq)) appears in Ref. |17|; it is monotonic and 
discontinuous everywhere. Thus, we can go directly from 
a qod orbit to a lower bound for the topological entropy. 
For example, the orbit R with code cn ~ 10011010 has 
r{R) = 1/3: any partially formed horseshoe which in- 
cludes an orbit of the braid type of R has entropy greater 
than ^1(^2/5) ~ 0.442. By taking rationals closer and 
closer to 1/3 we get better and better estimates (e.g., 
5/14 gives 0.481). We can also use the pseudo-Anosov 
orbits which are not qod to get estimates (e.g., the orbit 
s| mentioned above has h « 0.498), however in this in- 
stance we must compute the entropy from a train track 
and this can be a difficult calculation. 

In Table || we collect together some useful facts about 
braids in the horseshoe up to period eight. All the re- 
ducible orbits up to period nine only have finite order 
components. This is because the lowest period with a 
pseudo-Anosov orbit is five, and numbers less than or 
equal to nine do not have proper factors greater than 
or equal to five. As mentioned previously, the topologi- 
cal entropy for the pseudo-Anosov orbits which are not 
quasi-one-dimensional is calculated by finding the "train 
track" using the method of Bestvina and Handel H . This 
method results not only in a topological entropy, but also 



an explicit Markov partition (a 'topological model') and 
associated Perron- Frobenius matrix which is useful for lo- 
cating in phase space where the predicted periodic orbits 
are to be found. 



III. ROSSLER BRAID ANALYSIS 

A braid analysis of a low dimensional chaotic time 
series consists of four steps once an appropriate three- 
dimensional space is created [|l^: (i) the periodic orbits 
are extracted by the method of close recurrence [^,0 , 
(ii) the braid type of each periodic orbit is identified and 
the orbits are ordered by their two-dimensional forcing 
relationship [n,17| (iii) a subset of braids are selected 
which have maximal forcing and which force the orbits 
extracted in step (i), and (iv) if possible, an attempt is 
made to verify that some of the predicted orbits (not 
originally extracted in step (i)) are found in the system. 

In practice, steps (i) and (ii) are greatly simplified if 
the template or knot-holder organizing the flow can be 
identified using the procedure described by Mindlin and 
co-workers y,|l^,^,^ . Knowledge of the template helps 
in obtaining the symbolic names of the periodic orbits 
and in calculating the forcing relationship for the specific 
braids in that template. For instance, if the template is 
identified as a two-branch horseshoe knot holder (as are 
all the examples studied in this paper), then the theory 
of qod orbits of section II can be applied to simplify the 
analysis. 

Although template identification is very valuable, it is 
not essential for a braid analysis. Nor is the symbolic 
identification of the extracted orbits. In the worst case 
a braid analysis does require that the the braid conju- 
gacy class of each extracted periodic orbit is identified 
(see Elrifai and Morton |Q, or Jaquemard |Q for algo- 
rithms), and that the minimal Markov model (a 'train 
track' in the language of Thurston) can be constructed 
for each braid (see Bestvina and Handel H] , Los ||^ , and 
Franks and Misiurewicz |^ for algorithms). Algorithms 
exist for both of these steps, although the most computa- 
tionally efficient version of the braid conjugacy algorithm 
is probably not an effective solution beyond J5g. 

To illustrate the braid analysis we consider a chaotic 
attractor of the Rossler equations, 

X = ~{y + z) 

y ^ X + ey 

z = f + xz — iiz 

with e = 0.17, / = 0.4 and fi = 0.85. The Rossler equa- 
tions are integrated through 10^ cycles and the return 
map is examined at the half plane S = {{x,y,z) : x < 
& 2/ = 0}. As shown in Fig. 1, the template is easily 
identified as a horseshoe with zero global torsion. This 
template identification is verified by calculating the rela- 
tive rotation rates and linking numbers of the extracted 
periodic orbits as described by MindHn et. al. |15|,p9|,pl[ . 



To extract the (surrogate) periodic orbits by the 
method of close recurrence we first convert the return 
map from the sequence of values (a;„,z„) directly into 
the symbol sequence of O's and Vs. In this particular in- 
stance, since the map is close to one-dimensional, a good 
symbolic partition is obtained by examining the maxi- 
mum value of the next return map formed from the pro- 
jection on the x-coordinate — (— a;„, —x„+i) — at the 
surface of section. Orbits passing to the left of the maxi- 
mum are labeled zero, and those to the right are labeled 
one. Next we search this symbolic encoding for each and 
every periodic symbol string. Every time a periodic sym- 
bol string is found we calculate its recurrence and then 
save the instance of the orbit with the best recurrence. 
For instance, in searching for the period three orbit '100' 
we search the symbolic encoding of the return map for 
any instance of '100' and its cyclic permutations '010' 
and '001', and every time this symbol string is found, we 
next calculate its recurrence, which for this period three 
orbit is ey=o = (xn+s ~ Xn)'^ + {z„+3 - Zn)^, and then 
save the orbit with the minimum e. The advantage of 
this procedure of orbit extraction is that it is exhaustive. 
We search for every possible orbit up to a given period. 
In these studies we searched for all orbits between periods 
1 and 16. 

The resulting spectrum of periodic orbits up to period 
eight is shown in Table fS/\ This should be compared with 
the Table in Appendix A. The orbits which are present 
in (the full shift) Table of Appendix A, and not present 
in Table IV , are said to be pruned. Our goal is to predict 
as best as possible the pruned spectrum from the chaotic 
time series. 

Before we discuss the braid analysis, though, it is in- 
teresting to consider the number of orbits extracted as a 
function of the number of points in the return map. This 
is shown in Table ^ As expected, the number of orbits 
that can be extracted increases with the number of points 
in the return map. More importantly. Table M strongly 
suggests that using the method of close recurrence it is 
possible to obtain all the low period periodic orbits em- 
bedded within the strange attractor. For instance, after 
10* points are examined, we see that no new periodic 
orbits are found below period six. Similarly, after exam- 
ining 10^ points we believe we have found all orbits up to 
period eight. These results also caution us when working 
with small data sets — the extracted orbit spectrum is 
expected to miss orbits either because the orbit is pruned 
(it is not in the strange set) or because the sample of the 
strange set we are examining fails in providing a close 
enough coverage over the whole attractor. 

Using the results in section II, the finite-order and 
quasi-one-dimensional orbits in the extracted spectrum 
are easily identified from their orbit codes (see Table 
[V). Not unexpectedly, we see a sequence of quasi- 
one-dimensional orbits of increasing entropy (decreasing 
height) — the maximal qod orbit is the period 16 orbit 
1011011011011010 with entropy h w 0.480804. All orbits 
forced by this period 16 orbit up to period 8 are present. 



and none are missing. So this period 16 qod orbit already 
gives us a very good hyperbolic set with which to approx- 
imate to our (possibly nonhyperbolic) chaotic attractor. 
Can we do better? Doing better in this instance means 
identifying a pseudo-Anosov orbit which is not qod, but 
perhaps implies the maximal qod orbit found. Indeed, 
in this data set there is such an orbit, it is the period 8 
orbit 10010100 with entropy h w 0.498093. Again, the 
spectrum of orbits forced by this maximal pA orbit are 
consistent with the extracted spectrum which was exam- 
ined up to period 16. 

This data set is close to one-dimensional so kneading 
theory also does quite well for predicting the low period 
orbits. For instance, we could consider the period 3 orbit 
100, and this orbit (based on one-dimensional unimodal 
theory) also accurately predicts most of the extracted 
spectrum. However, this period 3 orbit is finite order, 
and we know from the results of Holmes and Whitley that 
there will be many (possibly high period) orbits which 
are forced by <i but not by <2 |3^ (in fact, we know 
that in 2-dimensions 100 forces only itself and a fixed 
point, and in 1-dimension it forces orbits of arbitrarily 
high period). Thus, although one-dimensional theory is 
a useful guide in this instance to the low period orbits, it 
can not be safely applied to make predictions about high 
period orbits. 



IV. BELOUSOV-ZHABOTINSKII REACTION 
BRAID ANALYSIS 

To further illustrate the braid analysis method we ob- 
tained data from the Bclousov-Zhabotinskii chemical re- 
action [Q . This is the same data set analyzed by Mindlin 
et. al. and consists of 65,000 equally spaced points which 
measure the time dependence of the bromide ion con- 
centration in the stirred chemical reactor. Following the 
techniques described by Mindlin et. al. we also embed- 
ded the scalar time series x{i), i = 1,2, . . . ,N in R"^ via 
a differential phase space embedding described by 



yi{i) ^x{i) 
2/2(i) = x{i) 
y^ii) = x{i) 



X*yi{i-1), A = 0.995 
x{i - 1) 



from which we reproduced the 3-dimensional attractor 
and return map shown in Figs. 5 and 6(a) of Ref. p5[ . 
The attractor is a zero global torsion horseshoe. There 
are approximately 125 points per cycle so the return map 
consists of about 520 points. 

Our technique for extracting (surrogate) periodic or- 
bits from this time series differs somewhat from that 
described in Ref. ||l^. We use the same procedure de- 
scribed for the Rossler data: first the return map data 
is converted into a symbol sequence of O's and I's de- 
pending on whether the orbit passes to the left or right 
of the maximum value of the return map, and second an 
exhaustive search is performed for all possible periodic 



orbits between periods 1 and 15. Again, this data set 
is almost one-dimensional and the simple symbolic pre- 
scription just described leads to a unique and consistent 
encoding of all the periodic orbits we are able to extract. 
Since the data set (at the return map) is small, we choose 
not to pick an arbitrary cut off for e, the close return cri- 
terion. Rather, we report the best e we are able to extract 
for a given periodic orbit (see Table VI). By including 
this additional piece of information we can make a more 
selective judgement about which close returns are, and 
are not, good surrogates for periodic orbits. In this way 
we are able to locate a few more periodic orbits than 

15|]. Also, one perhaps 



lis extraction method. 



were originally reported in Ref. 

surprising result comes out of t : 

It is not uncommon to find orbits of high period with 

very small recurrences. For example the period 13 orbit 

1011011101110 has a recurrence of e = 0.000712 which 

is significantly better than almost all of the lower period 

orbits. 

A list of the extracted orbits and their Thurston types 
is presented in Table VI . As expected, there is a sequence 



of quasi-one-dimensional orbits of increasing entropy, the 
largest of which is the period 16 orbit 1011011011011011 
with entropy h « 0.48084. In this particular instance, 
the period 16 qod orbit is in fact the maximal pseudo- 
Anosov orbit in the data set and it forces all the extracted 
orbits in this data set except for the finite order period 
three orbit 101. Indeed, a careful analysis of this data set 
suggests that the signal is subject to a small parametric 
drift which carries it between the strange attractor and 
a stable period 3 orbit (whose 1-dimensional entropy is 
/i« 0.4812) i7|. 



Table VII shows the number of predicted and extracted 



orbits as a function of period. The number of forced or- 
bits which are not in the extracted data set increases 
with the period. As with the Rossler data set, we be- 
lieve the forced orbits which are missing could actually 
be extracted from the data set if we were given a longer 
time series which provides a better coverage of the entire 
attractor. 



V. CONCLUSION 

In retrospect, we find it remarkable that such a small 
subset of periodic orbits (which are rather easy to get 
from experiments) contain so much topological and dy- 
namical information about a (low-dimensional) fiow. As 
previously demonstrated, a few low period orbits are suf- 
ficient to determine the template describing the stretch- 
ing and folding of the strange set MlqJsl] . The template 
provides an upper bound to the topological entropy and 
is, in a sense, a maximally (i.e., a full shift) hyperbolic 
set which we formally associate to a (possibly nonhy- 
perbolic) strange set p9| . In this paper we show how 
a sequence of periodic orbits (and their associated hy- 
perbolic sets) can be used to obtain a collection of finer 



and finer approximations to a strange set which is prob- 
ably not hyperbolic. Formally, we might say that the 
hyperbolic set associated with each pseudo-Anosov braid 
is embedded within the strange attractor we are trying 
to describe in the sense that the (possibly nonhyperbolic) 
strange set must contain at least all the orbits forced by 
the extracted pseudo-Anosov braid. We indirectly dis- 
cussed two measures of the goodness of this approxima- 
tion — the difference in the topological entropy, and the 
difference between the forced and extracted low period 
orbits. Using either measure, we have seen that it is pos- 
sible to select moderate period orbits (say < period 20) 
which provide good hyperbolic sets with which we can 
approximate a (possibly nonhyperbolic) strange set. 

The dynamical information derived from an orbit de- 
pends only on its braid type. As mentioned in section 
III, the braid type of an orbit can be determined without 
obtaining a good symbolic description of the orbits in the 
flow. We will illustrate these techniques in a future paper 
in which we consider a braid analysis of the bouncing ball 
system in a low-dissipation regime [ p8[ . Our strategy for 
handling the cases where a good symbolic description is 
not easy to obtain is to find a complete set of braid invari- 
ants on the small set of braids of interest. For instance, 
as mentioned in Table |, the exponent sum (simply the 
sum of the crossings) is a complete braid type invariant 
for horseshoe braids up to period 8. 

We would also like to point out that it is common for 
a collection of finite order braids to be pseudo-Anosov. 
This suggests an alternative strategy: instead finding a 
single orbit with maximal implications, it should also be 
possible to find a collection of orbits (possibly with very 
low period) which force all the observed orbits and pro- 
vides yet another hyperbolic approximating set. This is 
also the subject of future investigations. 
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APPENDIX A: APPENDIXES 

This Appendix shows the growth of periodic orbits in 
the unimodal map in Table VIIl and the unimodal or- 
dering in Table IX 
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FIG. 1. Orbits in the Rossler Attractor (e = 0.17, f = 0.4, 
jj, = 8.5): (a) Strange attractor showing surface of section, 

(b) the 1 and 10 orbits extracted from the chaotic time series, 

(c) the pretzel knot 1011010 extracted from the chaotic time 
series, (d) the pretzel knot shown as a braid. 



TABLE I. Low period horseshoe orbits. The exponent sum is a complete braid type invariant 
for aU horseshoe orbits up to period eight (see Ref. [17] for the exphcit conjugacy). The braid 
name notation is from Ref. [35]. Thurston braid types: finite order (fo), reducible (red), and 
pseudo-Anosov (pA), and quasi-one-dimensional (qod). The "red, fo" orbits are reducible Thurston 
types with finite order components. 



p 


cp 


■Kp 


Type 


es 


liP) 


h 


braid name 


symbolic name 


permutation 


Thurston type 


exponent sum 


height 


topological entropy 


Sl 


1 


(1) 


fo 





1/2 





S2 


10 


(12) 


fo 


1 


1/2 





S3 


w'i 


(123) 


fo 


2 


1/3 





/2x2 


1011 


(1324) 


red, fo 


5 


1/2 





S4 


lOO'i' 


(1234) 


fo 


3 


1/4 





si 


101 ly 


(13425) 


fo 


8 


2/5 





sl 


1001? 


(12435) 


pA, qod 


6 


1/3 


0.543535 


si 


1000? 


(12345) 


fo 


4 


1/5 





4 


10111? 


(143526) 


red, fo 


13 


1/2 





/3x2 


100101 


(135246) 


red, fo 


9 


1/3 





sl 


10011? 


(124536) 


red, fo 


9 


1/3 





sl 


10001? 


(123546) 


pA, qod 


7 


1/4 


0.632974 


si 


10000? 


(123456) 


fo 


5 


1/6 





4 


101111? 


(1453627) 


fo 


18 


3/7 





4 


101101? 


(1462537) 


pA, qod 


16 


2/5 


0.442138 


sl 


100101? 


(1362547) 


pA 


14 


1/3 


0.476818 


st 


100111? 


(1254637) 


pA 


14 


1/3 


0.476818 


S7 


100110? 


(1356247) 


fo 


12 


2/7 





s'r 


100010? 


(1246357) 


pA 


10 


1/4 


0.382245 


4 


100011? 


(1235647) 


pA 


10 


1/4 


0.382245 


4 


100001? 


(1234657) 


pA, qod 


8 


1/5 


0.666213 


4 


100000? 


(1234567) 


fo 


6 


1/7 





/22x2 


10111010 


(15472638) 


red, fo 


23 


1/2 





4 


1011111? 


(15463728) 


red, fo 


25 


1/2 





4 


101 1011? 


(14725638) 


fo 


21 


3/8 





4 


1001011? 


(13725648) 


pA 


19 


1/3 


0.346034 


4 


1001010? 


(13647258) 


pA 


17 


1/3 


0.498093 


si 


1001110? 


(13657248) 


pA 


17 


1/3 


0.498093 


st 


1001111? 


(12564738) 


pA 


19 


1/3 


0.346034 


sl 


1001101? 


(12573648) 


pA 


17 


1/3 


0.498093 


/2x4 


10001001 


(13572468) 


red, fo 


13 


1/4 





sl 


1000101? 


(12473658) 


pA 


15 


1/4 


0.568666 


st 


1000111? 


(12365748) 


pA 


15 


1/4 


0.568666 


4° 


1000110? 


(12467358) 


red, fo 


13 


1/4 





4^ 


1000010? 


(12357468) 


pA 


11 


1/5 


0.458911 


4^ 


1000011? 


(12346758) 


pA 


11 


1/5 


0.458911 


4' 


1000001? 


(12345768) 


pA, qod 


9 


1/6 


0.680477 


4^ 


1000000? 


(12345678) 


fo 


7 


1/8 






TABLE II. Topological entropy of quasi- one- dimensional 
(qod) orbits in a horseshoe up to period sixteen. 



Period 


Height 


Name 


Entropy 


5 


1/3 


imi'l 


0.543535 


6 


1/4 


lOOOli' 


0.632974 


7 


2/5 


101101^' 


0.442138 


7 


1/5 


100001? 


0.666213 


8 


1/6 


lOOOOOli' 


0.680477 


9 


3/7 


10111101?' 


0.397081 


9 


2/7 


10011001? 


0.604904 


9 


1/7 


10000001? 


0.687025 


10 


3/8 


lOUOUOl? 


0.473404 


10 


1/8 


100000001? 


0.690145 


11 


4/9 


1011111101? 


0.373716 


11 


2/9 


1000110001? 


0.65548 


11 


1/9 


1000000001? 


0.691662 


12 


3/10 


lOOUOUOOl? 


0.600332 


12 


1/10 


10000000001? 


0.692409 


13 


5/11 


lOlUlllUOl? 


0.361049 


13 


4/11 


lOUOUOUOl? 


0.479458 


13 


3/11 


1001 1001 1001? 


0.609 


13 


2/11 


100001100001? 


0.675816 


13 


1/11 


100000000001? 


0.692779 


14 


5/12 


lOllUOllUOl? 


0.41211 


14 


1/12 


1000000000001? 


0.692964 


15 


6/13 


lOUllllllUOl? 


0.354176 


15 


5/13 


10110111101101? 


0.467734 


15 


4/13 


10011011011001? 


0.599566 


15 


3/13 


10001 1001 10001? 


0.654699 


15 


2/13 


10000011000001? 


0.684869 


15 


1/13 


10000000000001? 


0.693056 


16 


5/14 


lOllOUOUOUOl? 


0.480804 


16 


3/14 


10001 10001 10001? 


0.656227 


16 


1/14 


100000000000001? 


0.693101 
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TABLE III. A sequence of quasi- one- dimensional (qod) or- 
bits with increasing forcing implications (entropy). Thus, ev- 
ery orbited forced by the period 7 qod orbit is also forced by 
the period 10 orbit, and so on. The notation "(*)" after an 
orbit indicates that one (but perhaps not both) of the sad- 
dle-node pair is forced. Both saddle-node partners will be 
forced by the next orbit in the qod sequence. 



Period 


Forced by (7) 
101101? 


Forced by (10) 
101 101 101? 


Forced by (13) 
lOUOUOUOl? 


1 








2 


10 






3 








4 


1011 






5 


1011?(*) 






6 


10111? 






7 


101111? 


101101? 




8 


10111010 
1011111? 


1011011?(*) 




9 


10111111? 
10111101? 
10110101? (*) 


101 101 11? 




10 


101111111? 
101111101? 
101110101? 
101101111?(*) 


lOUOlOlU 
101 101 101? 


lOUOUOl? 


11 


1011111111? 

1011111101? 

1011110111? 

1011110101? 

1011011111?(*) 

1011011101? (*) 


101 10101 11? 
1011010101? 


1011011011?(*) 
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TABLE IV. Periodic orbits extracted from time series data 
of the Rossler equations. All extracted orbits up to pe- 
riod 8 are shown. Between periods 9 and 16 only the ex- 
tracted quasi-one-dimensional (qod) are shown. The maximal 
pseudo-Anosov orbit (lOOlOlO'i) is not qod, but it forces all 
qod orbits with height greater than 1/3. Thurston Types: 



finite order (fo), reducible (red), and pseudo-Anosov (pA). 


Period 


Name 


Height 


Entropy 


Type 


1 


1 


1/2 





fo 


2 


10 


1/3 





fo 


3 


101 


1/3 





fo 


3 


100 


1/3 





fo 


4 


1011 


1/2 





red 


5 


10111 


2/5 





fo 


5 


10110 


2/5 





fo 


6 


101110 


1/2 





red 


6 


101111 


1/2 





red 


6 


100101 


1/3 





red 


7 


1011111 


3/7 





fo 


7 


1011110 


3/7 





fo 


7 


1011010 


2/5 


0.442138 


pA, qod 


7 


lOUOU 


2/5 


0.442138 


pA, qod 


7 


1001011 


1/3 


0.476818 


pA 


7 


1001010 


1/3 


0.476818 


pA 


8 


10111010 


1/2 





red 


8 


10111110 


1/2 





red 


8 


10111111 


1/2 





red 


8 


lOUOlU 


3/8 





fo 


8 


10110110 


3/8 





fo 


8 


10010110 


1/3 


0.346034 


pA 


8 


10010111 


1/3 


0.346034 


pA 


8 


10010101 


1/3 


0.498093 


pA 


8 


10010100 


1/3 


0.498093 


pA 


9 


101111010 


3/7 


0.397081 


pA, qod 


9 


lOllUOU 


3/7 


0.397081 


pA, qod 


10 


lOUOUOlO 


3/8 


0.473404 


pA, qod 


10 


lOUOUOU 


3/8 


0.473404 


pA, qod 


11 


10111111010 


4/9 


0.373716 


pA, qod 


13 


lOUOUOUOlO 


4/11 


0.479450 


pA, qod 


13 


lOllOUOUOU 


4/11 


0.479450 


pA, qod 


15 


lOUllllllllOU 


6/13 


0.354176 


pA, qod 


15 


lOllOllUOUOlO 


5/13 


0.467734 


pA, qod 


16 


lOllOUOUOUOlO 


5/14 


0.480804 


pA, qod 
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TABLE V. Number of periodic orbits extracted (e = 0.005) 
as a function of the number of points in the return map for 
time series data from the Rossler equation. The data suggests 
that after 100,000 points, all the periodic orbits in the time 
series up to period eight are extracted. 



Period 


100 


500 


1000 


10,000 


50,000 


100,000 


1 





1 


1 


1 


1 


1 


2 





1 


1 


1 


1 


1 


3 








2 


2 


2 


2 


4 








1 


1 


1 


1 


5 











2 


2 


2 


6 


1 


1 


1 


3 


3 


3 


7 


1 


1 


2 


6 


6 


6 


8 


1 


2 


2 


8 


9 


9 


9 





2 


4 


9 


14 


14 


10 





1 


4 


11 


17 


19 
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TABLE VI. Periodic orbits extracted from the Be- 
lousov-Zhabotinskii Reaction time-series up to period 15. All 
orbits with a best (normalized) recurrence of less than 0.1 are 
shown. The period 16 orbits are from Ref. [15]. 



Period 


Name 


Recurrence 


Type 




1 


1 


0.016782 


fo 


2 


10 


0.002615 


fo 


3 


101 


0.000128 


fo 


4 


1011 


0.002648 


fo 


5 


10111 


0.002962 


fo 




5 


10110 


0.013668 


fo 




6 


101110 


0.006449 


fo 




6 


101111 


0.029014 


fo 




7 


1011110 


0.005088 


fo 




7 


1011010 


0.041837 


pA, qod, h - 


0.442138 


7 


1011011 


0.010585 


pA, qod, h - 


0.442138 


8 


10111010 


0.014287 


fo 




8 


10110111 


0.017370 


fo 




8 


10110110 


0.070293 


fo 




9 


101111010 


0.001720 


pA, qod, h - 


0.397081 


9 


101101011 


0.061843 


pA 




9 


101101010 


0.020884 


pA 




9 


101101110 


0.018380 


pA 




9 


101101111 


0.050094 


pA 




10 


1011101010 


0.037910 


fo 




10 


1011101011 


0.007237 


fo 




10 


1011111010 


0.035966 


fo 




10 


1011011010 


0.003440 


pA, qod, h - 


0.473404 


10 


1011011011 


0.008971 


pA, qod, h - 


0.473404 


11 


10110101010 


0.038044 






11 


10110111010 


0.009119 






11 


10110110111 


0.008447 






11 


10110110110 


0.085593 


fo 




12 


101110101011 


0.031644 






12 


101110101010 


0.013333 






12 


101101010110 


0.029674 






12 


101101111010 


0.055765 






12 


101101101011 


0.032222 






12 


101101101010 


0.016560 






12 


101101101110 


0.048563 






12 


101101101111 


0.064420 






13 


1011110101010 


0.052853 






13 


1011010111010 


0.022397 


pA 




13 


1011010101110 


0.004291 






13 


1011011101110 


0.000712 






13 


1011011101010 


0.011199 


pA 




13 


1011011111010 


0.046255 


pA 




13 


1011011010110 


0.097165 


pA, qod, h - 


0.479450 


13 


1011011011010 


0.013063 


pA, qod, h - 


0.479450 


13 


1011011011011 


0.004663 






14 


10110101110110 


0.056515 






14 


10110101010110 


0.009366 






14 


10110111010110 


0.047037 






14 


10110110111010 


0.060024 






14 


10110110110110 


0.004920 


fo 




15 


101101011101110 


0.034250 






15 


101101110101010 


0.015669 






15 


101101101011010 


0.010945 






15 


101101101010110 


0.076061 







14 



15 


101101101110110 


0.049261 


pA 




15 


101101101111010 


0.007018 






15 


101101101101010 


0.049573 






15 


101101101101110 


0.010284 






15 


101101101101111 


0.078043 






16 


1011011011101010 


na 


pA 




16 


1011011011111010 


na 


pA 




16 


1011011011011010 


na 


pA, qod, h - 


= 0.480804 


16 


1011011011011011 


na 


pA, qod, h - 


= 0.480804 



TABLE VII. Number of periodic orbits extracted and pre- 
dicted for time series data from the Belousov-Zhabotinskii 
reaction (fa 500 points in the return map). As illustrated 
with the Rossler data, we expect that the time series is far 
to short to be able to extract all the predicted orbits of any 
except the lowest periods. 



Period 


# predicted 


# found 


^ missing 


^ marginal 


1 


1 


1 






2 


1 


1 






3 


1 


1 






4 


1 


1 






5 


2 


2 






6 


2 


2 






7 


4 


3 


1 




8 


5 


2 


2 


1 


9 


8 


4 


2 


2 


10 


11 


5 


6 





15 



TABLE VIII. The number of periodic orbits in a horseshoe 
of least period n is given given by (^^ 2'^(f>{n/d))/n where d 
ranges over all the divisors of n and (j> is the Mobius func- 
tion defined by <^(m) = if m has a repeated prime factor, 
otherwise 0(m) = (_ !)("'""''<='■ "f p"™« /"'=t°'-» "f ™). 





Number of Periodic Orbits 


Period 


of Least Period 


1 


2 


2 


1 


3 


2 


4 


3 


5 


6 


6 


9 


7 


18 


8 


30 


9 


56 


10 


99 


11 


186 


12 


335 


13 


630 


14 


1161 


15 


2182 


16 


4080 


17 


7710 


18 


14532 


19 


27594 


20 


52377 


21 


99858 


22 


190557 


23 


364722 


24 


698870 


25 


1342176 


26 


2580795 


27 


4971008 


28 


9586395 


29 


18512790 


30 


35790267 
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TABLE IX. One-dimensional unimodal ordering up to pe- 
riod 10. Orbits are sorted first by period and second by uni- 
modal ordering. Orbits agreeing in all but the last digit (0,1) 
arc braid equivalent saddle-node pairs. 



Period 


Name 


1 





1 


1 


2 


10 


3 


101 


3 


100 


4 


1011 


4 


1001 


4 


1000 


5 


10111 


5 


10110 


5 


10010 


5 


10011 


5 


10001 


5 


10000 


6 


101110 


6 


101111 


6 


100101 


6 


100111 


6 


100110 


6 


100010 


6 


100011 


6 


100001 


6 


100000 


7 


1011111 


7 


1011110 


7 


1011010 


7 


1011011 


7 


1001011 


7 


1001010 


7 


1001110 


7 


1001111 


7 


1001101 


7 


1001100 


7 


1000100 


7 


1000101 


7 


1000111 


7 


1000110 


7 


1000010 


7 


1000011 


7 


1000001 


7 


1000000 


8 


10111010 


8 


10111110 


8 


10111111 


8 


10110111 


8 


10110110 


8 


10010110 


8 


10010111 


8 


10010101 


8 


10010100 


8 


10011100 


8 


10011101 


8 


10011111 


8 


10011110 


8 


10011010 



8 


10011011 


8 


10001001 


8 


10001011 


8 


10001010 


8 


10001110 


8 


10001111 


8 


10001101 


8 


10001100 


8 


10000100 


8 


10000101 


8 


10000111 


8 


10000110 


8 


10000010 


8 


10000011 


8 


10000001 


8 


10000000 


9 


lOllUlU 


9 


101111110 


9 


101111010 


9 


lOllUOU 


9 


101101011 


9 


101101010 


9 


lOUOlUO 


9 


101101111 


9 


100101100 


9 


100101101 


9 


100101111 


9 


100101110 


9 


100101010 


9 


100101011 


9 


lOOlUOU 


9 


100111010 


9 


100111110 


9 


100111111 


9 


100111101 


9 


100111100 


9 


100110100 


9 


100110101 


9 


lOOUOlU 


9 


lOOUOUO 


9 


100110010 


9 


lOOUOOU 


9 


100010011 


9 


100010010 


9 


100010110 


9 


100010111 


9 


100010101 


9 


100010100 


9 


100011100 


9 


100011101 


9 


100011111 


9 


100011110 


9 


100011010 


9 


lOOOUOU 


9 


100011001 


9 


100011000 


9 


100001000 


9 


100001001 


9 


100001011 


9 


100001010 



17 



9 


100001110 


9 


100001111 


9 


100001101 


9 


100001100 


9 


100000100 


9 


100000101 


9 


100000111 


9 


100000110 


9 


100000010 


9 


100000011 


9 


100000001 


9 


100000000 


10 


1011101010 


10 


lOlUOlOU 


10 


lOlllUOU 


10 


1011111010 


10 


lOllUlUO 


10 


lOlUlllU 


10 


lOUOlOlU 


10 


lOUOUlU 


10 


1011011110 


10 


1011011010 


10 


lOUOUOU 


10 


lOOlOUOU 


10 


1001011010 


10 


1001011110 


10 


1001011111 


10 


1001011101 


10 


1001011100 


10 


1001010100 


10 


1001010101 


10 


1001010111 


10 


1001010110 


10 


1001110010 


10 


lOOlUOUO 


10 


lOOlUOlU 


10 


1001110101 


10 


1001110100 


10 


1001111100 


10 


1001111101 


10 


lOOllUlU 


10 


1001111110 


10 


1001111010 


10 


1001111011 


10 


1001101011 


10 


1001101010 


10 


lOOUOlUO 


10 


1001101111 


10 


1001101101 


10 


lOOUOUOO 


10 


1001100100 


10 


1001100101 


10 


lOOUOOlU 


10 


1001100110 


10 


1000100110 


10 


1000100111 


10 


1000100101 


10 


1000100100 


10 


1000101100 


10 


1000101101 



10 


1000101111 


10 


1000101110 


10 


1000101010 


10 


1000101011 


10 


1000101001 


10 


1000101000 


10 


1000111000 


10 


1000111001 


10 


lOOOlUOU 


10 


1000111010 


10 


1000111110 


10 


1000111111 


10 


1000111101 


10 


1000111100 


10 


1000110100 


10 


1000110101 


10 


lOOOUOlU 


10 


lOOOUOUO 


10 


1000110010 


10 


lOOOUOOU 


10 


1000010001 


10 


1000010011 


10 


1000010010 


10 


1000010110 


10 


1000010111 


10 


1000010101 


10 


1000010100 


10 


1000011100 


10 


1000011101 


10 


1000011111 


10 


1000011110 


10 


1000011010 


10 


lOOOOUOU 


10 


1000011001 


10 


1000011000 


10 


1000001000 


10 


1000001001 


10 


1000001011 


10 


1000001010 


10 


1000001110 


10 


1000001111 


10 


1000001101 


10 


1000001100 


10 


1000000100 


10 


1000000101 


10 


1000000111 


10 


1000000110 


10 


1000000010 


10 


1000000011 


10 


1000000001 


10 


1000000000 



18 



This figure "bra_fig_l.gif" is available in "gif ' format from: 



http://arXiv.org/ps/chao-dyn/9411008v2 



